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Abstract 

In this paper, we propose a method based on GMM (the generalized method 
of moments) to estimate the parameters of stable distributions with < a < 2. 
We don't assume symmetry for stable distributions. 
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1 Introduction 

The class of stable distributions is an important probability distribution class, 
which can be viewed as the sum of a large number of independent and identi- 
cally distributed random variables with very small individual effects. The stable 
non-Gaussian distributions has the property that the second moment is infinite. 
Plenty of evidence suggested that some important economic variables such as stock 
price changes, interest rate changes, and price expectations etc. can be better de- 
scribed by stable distributions, in most cases, by stable non-Gaussian distributions, 
see Mandelbrot (1963a,b), Fama (1965) etc.. It was suggested by Bartels (1977), 
Koenker and Bassett (1978) that the distribution of the regression disturbance may 
also belong to the class of stable distributions 
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One dimensional stable distributions with parameter a = 1 and a = 2 are 
Cauchy distribution and normal distribution respectively. Both are well studied. 

For stable distributions with < a < 2, since there are no closed form ex- 
pressions for the density function in most cases, estimation of parameters encoun- 
ters difficulties. Some important works in this area were done by Fama and Roll 
(1968,1971) and others. Stable distributions were also used in tests for nomality by 
Bera and McKenzie (1986). 

Previous methods on estimation of parameters of stable distributions include the 
fractile method of Famma and Roll (1968,1971), the improved version of McCul- 
loch (1986), the approximate maximum likelihood theory developed by DuMouchel 
(1973a, b, 1975), and the iterative regression method of Koutrouvelis (1980,1981) etc.. 
Akgiray and Lamoureux (1989) made a comparative study of the fractile method 
and iterative regression method. 

In this paper we propose a method based on GMM (the generalized method of 
moments) to estimate simultaneously all the parameters of stable distributions with 
1 < a < 2 and < a < 1. We don't assume the symmetry of stable distribution 
here. 

Stable distributions with < a < 2 and a^l may be defined by the character- 
istic function 

ITtt 

jj, = exp[— c|z| a (l — i/3tan— sgnz) + irz], (1.1) 

with c > 0, (3 € [—1, 1] and r € R, where tan is tangent function and sgn is sign 
function. (3 = corresponds to symmetric stable distributions. (3^0 corresponds to 
nonsymmetric stable distributions. Two special cases (3 = 1 and (3 = — 1 correspond 
to one sided stable distributions. Here we focus on the cases that — 1 < (3 < 1. r 
is the drift term. For detailed account of stable distributions, see Sato (1999) and 
Zolotarev (1983). 

We assume the random variable x has the stable distribution with parameters 
a, (3, t, and c. 

In order to present infinite and finite series expressions of the density function 
of the stable distribution in Sato (1999), we introduce a new parameter c which 
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satisfies 



2 — a 

c = cos(7r/3 — - — )c a , if 1 < a < 2; 
c = cos(^— )c a , if < a < 1. 



Thus estimation of c is equivalent to estimation of c. Consider the new variable 



X — T 
U = . 



u has cos(7r /3^p) if 1 < a < 2 and cos(^y^) if < a < 1 as a parameter which 
appears in the position of c in the characteristic function expression (1.1). 

By Sato (1999), we have the following infinite and finite series expressions of the 
density function associated with variable u. 

1. When 1 < a < 2, the convergent series expression for the density function is 

Pa ( u ) = - ^(_ 1 )n-i r ( TO /« + l) ( sin ( 7m ( 1 _ pl^Lj/Vfiu 1 ^ 1 , for iieR, 

n=l 

(1.2) 

and for u — > oo, the finite series expression is: 

p a ( u ) = I y ( _ 1) n-i r (^ + 1 ) (sin(7rn(1 _ ^i^)/ 2a )) u -na-l 

^ y ' vr ^ V ; n\ K K K a " " (13) 

rz=l v y 

+ 0(u-^+ 1 )"- 1 ). 

For -u — > — oo, since the dual distribution u = —u has the density p a (u) = 
Pa(—u), the finite series expression is: 

p>) = ^(-ir 1 ^^ (L4) 

+ 0(\u\-^ a - 1 ). 

2. When < a < 1, for u > 0, the convergent series expression for the density 
function is 

p a (u) = - jr(-l)^ T{na ± 1} (sin(7rn(l + 0)/2a)) U -™-\ (1.5) 
7r n! 

n=l 
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and for u < 0, the convergent series expression for the density function is 

p a (u) = - f;(-l)"- i r( " a + 1) (sin(7m(l - 0)/2a))\u\-™-\ (1.6) 
7r z — ' n! 

n=l 

which is obtained by duality argument similar to the above. 
For u — > 0, the finite series expression is: 

p a (u) = - V(-ir- i r(re/a , + 1) (sin(vrn(l + / 9)/2))«"- 1 + 0(u"). (1.7) 
7T ' n! 

n=l 

In the following two sections, we discuss a method based on GMM to estimate 
parameters a, (3, r, and c. For convenience, we denote the vector (a, (3, r, c) by 6 and 
the vector (ao, /?o> r o> Co) by #o- in section 2, we look at the case when 1 < a < 2. 
In section 3, we study the case when < a < 1. 



2 1 < a < 2 



Let us consider the following moment functions: 



fc+i 



/*(y;0) = U 



»2(1+ 



2(fc+l) 



rf 2 (^)^: 



if y>i?i; 



if i? 2 < y < Ri; 



if y <R 2 , 



(2.1) 



for fc 



1,2, ...,m with m > 4, where the exponent 2(1 + 2 (k+i) ) means nrs t take 



square and then take exponent 1 + 



2(fe+l) • 



i?i and i?2 are constants such that 



R2 < t < R\ with ' Rl g ~' r and — ^ not large. The reason for this requirement is 
that the convergent formula (1.2) of the density function is appropriate for small 
u and the finite series expressions (1.3) and (1.4) are appropriate for u is large. 
This phenomenon was pointed out by Famma and Roll (1968) for symmetric stable 
distribution. When we solve the optimization problem resulted by GMM, if the 
solution for r is very close to R\ or R2, it means that we may not get the globally 
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optimization yet. We have to make R\ bigger or R 2 bigger, and do optimization 
again with the solution from previous step as the initial point. About d\ and d 2 in 
the above expression, they are constants satisfying 

d-y = {{Ri - t)/c) 2 and d 2 = ((r - R 2 )/c) 2 . 

We can see that fk(y,9), k = 1, ...,m, are continuous functions. 

Let x be the random variable of stable distribution with parameters 9. Then 

E(f k (x;9)) = [°° f k (y-9)p a (K 1 )ldy, (2.2) 

J — oo c c 

where p a is the density function in (1.1), (1.2) and (1.3). 

Let xi, X2, x n be independent stable distributions with the same parameter 
vector 9q. Since E\fk(xf, 9)\ < 00, i = 1, ...,n, fe = 1, ...,m, by strong law of large 
numbers, we have 



n 

as n — > 00. 



E(fk(x;9o)) almost surely, fe = 1, ...,m, 



Denote the sample moment fki^tfi)) ^ E n fk(x;9), k = 1, ...,m. 
We can see that E n fk(x;9)) and Efk(x;9) have first derivatives (vector) with 
respect to = (a,(3,T, c). Define 

V = ci[£ n /i(x; (9 ) " ^))] 2 + ■ • • + c m [i? n / m (x; 0„) - S/ m (x; ^)] 2 , (2.3) 

where ci, ...,c m are positive constants. Then by GMM, the minimization solution 9 
of V converges in distribution to #o- 

Since we don't have closed form expression for the density of stable distributions, 
we need approximate finite series expressions for Efk(x;9)), fe = 1, ...,m. When 
we minimize the function V with respect to the parameter vector 9, we use the 
approximate finite series expressions for Efk(x; 9) instead of Efk(x; 9), k = 1, m 
themselves. 

In the following, we will show the existence of approximate finite series expres- 
sions for Efk(x; 9)), fe = 1, m. 
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For positive integer N > 0, define Pi,n, P2,N and P3,n as follows: 



p 1)N {u) = - ^(-lf - 1 ^^±H( s in(7rn(l - p^—^)/2a)) U - na -^ for u > 0, 

n=l 

(2.4) 

N 

P2,n(u) = -V(-l) n - i r(rea , + 1) (sin(vrn(l+/3^)/2«))H-""- 1 , for u < 0, 
7r ^— n! a 

n=l 

(2.5) 

P3,n(u) = - V(-l)"- i r(n/a + 1) ( S m(vrn(l - /J^)^))^" 1 , for M eK, 

n=l 

(2.6) 

withp^jv, P2,N andp3 5 Ar corresponding to (1.3), (1.4) and (1.2) respectively. 
By (1.3), we have 

f°° 1 
/ fk(v,0)\Pa((y ~r)/c) -pi, N ({y -T)/c)\-dy 
JRt c 

r°° i 

< Mi / u—u-^+^^du 

c 

( where M\ is a constant ) 

poo . 

< M 1 / u-^ a ~ 1+ ^Tdu 

c 

< 0((( J R 1 -r)/5)' (iV+1)a+ ^T), 

which goes to as N — > oo, when (i?i — r)/c > 1. 
Similarly, by (1.4), we have 



(2.7) 



/ fk{y,0)\Pa((y - t)/c) -P2,n((v - r)/c)\-dy 

J —oo c 

< 0(({T-R 2 )/5)- (N+1)a+ ^), 



(2.8) 



which goes to as N — > oo when (r — R 2 )/c > 1. 
By (1.2), we have 

Pi,n({v ~ t)/c) -> p a ((y - t)/c) pointwisely for y G [-R 2 ,-Ri], as iV -> oo, 

and the convergence is uniform when both ((i?i — r)/c) and ((r — -R 2 )/c) are less 
than 1. 
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In fact, we can show that for given 9 = (a, (3, r, c), given R\ and R±, without the 
assumption that ((R\ — r)/c) and ((r — R2)/c) are less than 1, 

Pi,n((v-t)/c) ^p a ({y-r)/c) uniformly for all y £ [R2, Ri], as N -> 00 (2.9) 

To show this statement, we need the following lemma: 
Lemma 2.1 /or given M > 0, L > 1, there exists a positive integer no, which 
depends on L, M and 1 < a < 2, such that when n > no, 

r(m/a , + 1) M"-' < L~ n . (2.10) 

Proof. We know that 

T(n/a + l) <C ([n/a] + l)\, (2.11) 

where [n/a] is the greatest integer less than n/a and Co is a positive constant. 
Since 

n' 

> 1, for m = 0,1, ...,n, (2.12) 



m!(n — m)! 
letting m = [n/a] + 1, then 

([n/a] + 1)! 



< 7 r / i iTT ' ( 2 - 13 ) 



n! (n — [n/a] — 1)1 

which means 

r(n/a , + 1) M- 1 < C ([n/a] , + 1)! M- 1 < C 07 r ^ , M" -1 . (2.14) 

n! n! (n — [n/aj — 1)! 

By Stirling's formula, 

(n-[n/a]-l)\ ~ (n- [n/a] - i)(n-[n/a]-i) e -{n-[n/a]-i) ^{/n~^ja] - 1), (2.15) 
as n — > 00. 

Thus, in order to show the lemma, it is sufficient to show that there exists no > 
such that for any n > no, 

C M n ~ l L n < (n - [n/a] - i)(n-[«/«]-i) e -("-[«/«]-i)^( n _ [„/ a ] _ 1). ( 2 .16) 
It is clear that 



log{(n - [n/a] - i)(n-[n/a]-i) e -(n-[n/ a ]-i) - [n/a] - 1)} 

1 

2 



(n - [n/a] - 1) log(n - [n/a] - 1) - (n - [n/a] - 1) + - log(27r) (2.17) 



+ ilog(n - [n/a] - 1), 



where (n — [n/a] — 1) log(n — [n/a] — 1) becomes the dominant term when n is large 
enough. 

We can see that when n is large enough, 



log Co + (n - 1) log M + n log L < (n - [n/a] - 1) log(n - [n/a] - 1), (2.18) 
which means that there exists no > such that for any n > no, (2.16) holds. 



Let M = min((i?i — r)/c, (r — R2)/c). Substituting for u in (1.2), applying 
Lemma 2.1 and using the following fact, 



we have (2.9) holds. 

One thing we need to pay attention is that although we have uniform convergence 
in (2.9), from the proof of Lemma 1, we can see that the convergence speed depends 
on (i?i — t)/c, (r — R2)jc and a. When a is close to 1 or either R\ or R2 is very 
large, the convergence speed will be slow. Since we cannot change the parameter a, 
t and c, the only way to obtain good convergence speed is to keep R\ and R2 close 
to r. 

Clearly (2.9) implies that 



Therefore we approved the lemma. 



□ 





(2.19) 



as N 



00. 



Define 




(2.20) 



for k = 1, m. 
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Combining (2.7, (2.8) and (2.19), we have 

T kll , N + T k ^ N + T kAN -> £7(/ fc (x; 0)), (2.21) 

for k = 1, ...,m, as iV goes to oo. This means the sum of T kj i^, T k ^,Ni and T k ^,N 
is the approximate finite series expression we wanted. 

Next we give the expressions of T k ^^, T ky2 ,N, and T k ^,N without integrals inside. 
It is easy to see 



1. 



T k ,l,N 

f°° 1 
/ fk(y;0)pi,N{(y -r)/c)-dy 

JRi c 

N 



(^ E( _ ir -,r( ! ^+i) (sinM1 _ (3 2-« )/2a)) 



(2.22) 



n=l 

"CO 



/■°° _i . 

• / u*+i?j na L du, 



where 



°° i 1 



roc 

= T u- na - 1+ —idu (2.23) 

c 



2. 



T k) 2,N 
ri?2 



= / fk{y,6)P2,N{(y -r)/c)-zdy 

J — oo ^ 

_ (fr-ft)/°Tf (■p.-.nna + l) + ^ )/2a)) (2-24) 

7T ^— ' n! a 

n=l 



i? 2 -T 

• / C \u\^\u\~ nOL ~ x du, 

J —CO 
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where 



3. 



where 



with 



and 



I C \u\T^\u\- na - l du 

J — ( 



oo 
oo 



u - («+!)+ TTT du 
t-R 2 

C 

1 ((r-R 2 )/c)- na+ ^. 



Tk,3,N 

/ fk(v,0)P3,N((y -r)/c)-dy 



iRi 

N 



1 i\n-i r ( n / a + 1) / • / /i a 2_a \/o\\ 

-Z^( _1 qsin7rnl-/3 )/2)) 

7T n! a 

n=l 



u {l+ ^) ] u n - l du, 



/" c ocih ,-i ) r , i , if n is odd: 

/ u A + ^+^'u n - 1 du = { 

J R K T I A\ — A2, if n is even, 



A l = I u 2(1+ 5^Fij) u »- 1 d u 



«1- 



/o 



1 /Y „\ /^+ 2 (H-2(fcW) 



ra + 2 (! + 2(iTi)) 



^-((^i-r)/c) 



^0 

— — ^-((r-^/c)^ 1 ^'. 
n + 2 ( X + 2(1+1) ) 
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Thus we obtain the expressions of Tfc^jv, T k> 2 t N, and Ifc^jv without integrals in- 
side. By (2.21), we can substitute the sum of T k; ± t N, Tk,2,N and T^n for E(f k (x; 9)) 
in the expression of F of (2.3) and then do optimization to get the estimate 9 for 
9 - 



3 < a < 1 



In order to ensure integrability, we consider the following moment functions, 
which are different from those when 1 < a < 2: 

-I" 

c 



fk(y;9) 



rfi(^)^ 1 , Hy>Ri; 

(^Ifd-OT)), ifi? 2 <y<i? i; (3.1) 



^(I^EJIS+IJ, ify<i? 2 , 

for fc = 1,2, ...,m with m > 4, where the exponent 2(1 — 2 (k+i) ) means first take 
square and then take exponent 1 — 2 (k+i) • ^ an< ^ ^ 2 are cons tants such that 
i?2 < t < i?i with Rl r T and r ~ g R2 are small. The reason for this requirement 
is that when u is small, we have the finite series expression (1.7) for the density 
function. d\ and d 2 in the above expression are constants satisfying 

di = {{Ri - t)/c) 2 and d 2 = ((r - i? 2 )/c) 2 . 

It is clear that f k (y,9), k = 1, ...,m, are continuous functions. 

Let x be the random variable of stable distribution with parameters 6*. Then 

E(f k (x-9)) = f°° f k ( y -9) Pa (K^)ldy, (3.2) 

7 — oo ^ ^ 

where p a is the density function in (1.5), (1.6) and (1.7). 

Similar to the argument in section 2, for independent stable distributions xi, x 2 , x n 
with the same parameter vector 9q, denote the sample moment — Ih^n?)} by 
E n fk(x;9), k = l, ...,m. Since E\f k (xf,9)\ < oo, i = 1, ...,n, fe = 1, ...,m, by strong 
law of large numbers, we have 

E n fk(x;9 ) -> E(f k (x;9 )) almost surely, fe = 1, ...,m, 
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as n — > oo. 

Both E n fk(x;6)) and Efk(x;0) have first derivatives (vector) with respect to 
= (a, (3, t, c). We define 

V = c^EnUx- Oo) - Efi(x; 0))} 2 + ■■■ + c m [E n f m (x; ) - Ef m {x; 0)} 2 , (3.3) 

where c\, ...,c m are positive constants. Then by GMM, the minimization solution 9 
of V converges in distribution to Oo- 

Because of no closed form expressions for the density of stable distributions, we 
substitute its approximate finite series expression (if they exist) for Efk(x; 0)), k = 
1, ...,m in the expression of V of (3.3). 

Using the similar argument as we did in section 2, it is not hard to show the 
existence of approximate finite series expression for Efk(x;0)), k = l,...,m. For 
N > 0, we denote it by the sum of Tk,i,N, T^^n, and T^n, where 



f°° - 1 
Tk,i,N= / fk{y,9)pi,N({y-T)/c)-dy, 
JRi c 

rR 2 _ i 

Tk, 2 ,N= / fk(y,0)P2,N((y-T)/c)-dy, (3.4) 

f R1 - 1 

T k ,3,N = / fk{y;0)P3,N({y -r)/c)-dy, 

JR->. c 



IR 2 

for k = 1, ...,m, with p ljN , p 2 ,N and p 3tN being: 
N 



Pl,N 



(«) = 1 ^(-i)"~ 1 £i!^±ll ( s i ri ( 7 rn(l + (3) /2a))u- na ' 1 , for u > 0, (3.5) 

n=l 

N 

p 2iN (u) = - y(-l) n - i r(rea , + 1) (sin(7rn(l - /3)/2a))|uj- na - 1 , for u < 0, (3.6) 
7r z — ' n! 

n=l 

P3,7v(n) = - V(-l)"- 1 rWa , + 1} (sin(7rn(l + 0)/2))u n -\ for u G R, (3.7) 



7T * — * n! 

n=l 



We can see that pi, jv, P2,at and P3 5 at corresponding to (1.5), (1.6) and (1.7) respec- 
tively. 

Next we give the expressions of T^^at, T^n, and Tfe^Ar without integrals inside. 
It is easy to see 
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1. 

Tk,l,N 

f°° - 1 
= / fk(v,0)Pi,N((v-T)/c)-dy 



where 



C 2 7T 7T * — ' ' n! 
ra=l 



oo 



ufc+in na 1 du 

C 

OO j 

_ u _na_1 ~fe+T(in (3.9) 



2. 



2JV 



R 2 i 

fk(y;0)p 2 ,N((y - r)/c)-dy 



— oo 



«r - ft)/* f ( _ ir .r(«. + i) (rio(m p _ W2a)) (3.io) 



7r ' — * n! 

n=l 

/ |n|fe+T|n|~ Tta_1 (iu, 



where 



f c W^WV^^du 
J — < 



-oo 

POO 



/>oo . 

/ u - {n+1) - — du (3.11) 



L_ ((r _ j R 2 ) /a) — sk. 
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3. 



where 



with 



and 



Tk,3,N 

/ fk(y;0)p3,N({y -T)/c)-d y 

JRi c 

l E(-ir i r(n/ ;, + 1) (sin (7 rn(l + fl/2)) (iU2) 



n=l 



= Mi + A 2> if n is odd; 



I >li — ^2, if n is even, 



, 2 ( 1 ~2TFFiT)„n-l / 



JO 

V ^-^-WW^du (3-14) 



I 



ra + 2( 1 "2(ITi)) 



((fli - T )/c) n+2(1 "™y ) 



-«2 

, 2 ( 1— 2(fcTT) )<j/ n— -*v 



M »-'+*-^l> dM (3.15) 



-«2 



- j— ((r - R 2 )/c) n+2{1 ~^ 



Now we have the expressions of Tjfc^jv, 7fc,2,iv j an d T\ j3) jv without integrals inside. 
We can substitute the sum of X^i^/v, Tk y2 ,N and Tfc^jv for E(fj c (x;6)), k = 1, ...,m 
in the expression of V" of (3.3) and then do optimization to get the estimate for 
O . 
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